Pairing symmetry signatures of T\ in superconducting ferromagnets 
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We study the nuclear relaxation rate X/Tj as a function of temperature for a superconducting- 
ferromagnetic coexistent system using a p-wave triplet model for the superconducting pairing sym- 
metry. This calculation is contrasted with a singlet s-wave one done previously, and we see for the 
s-wave case that there is a Hebel-Slichter peak, albeit reduced due to the magnetization, and no 
peak for the p-wave case. We then compare these results to a nuclear relaxation rate experiment on 
UGe2 to determine the possible pairing symmetry signatures in that material. It is seen that the 
experimental data is inconclusive to rule out the possibility of s-wave pairing in UGei. 
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In the BCS theory of superconductivity, the conduc- 
tion electrons in a metal cannot be both ferromagnet- 
ically ordered and superconducting. Superconductors 
may expel small magnetic fields passing through them, 
but strong magnetic fields destroy the superconductivity 
(SC). Even small amounts of magnetic impurities are usu- 
ally enough to eliminate SC. Much work has been done 
both theoretically and experimentally to understand this 
interplay, and to explore the possibility of coexistence be- 
tween these two ordered states. 

The archetypal material that does display this coexis- 
tence is UGe2, and it is generally believed that the su- 
perconductivity is governed by triplet pairing due to the 
presence of the internal ferromagnetism 2 . However, there 
have been some convincing calculations that instead con- 
sider s-wave superconductivity to be favored in the pres- 
ence of this magnetization. For instance, Abrikosov 3 
and Suhli demonstrate that XJGe% should have s-wave 
Cooper pairing which is mediated by localized ferromag- 
netic spins. Also, Blagoev et al. 5 studied a weak fer- 
romagnetic Fermi liquid and showed that s-wave super- 
conductivity is possible and subsequently energetically 
favored in the ferromagnetic state. 

The contributions of all of these theories and calcu- 
lations have naturally led to a healthy debate on this 
specific property of this fascinating material, and thus 
the nature of the pairing is still an open question™. 

Motivated by this debate and the experimental data 
that should hold the answer, the present authors 
have recently studied various properties of a mean-field 
model of coexistent itinerant ferromagnetism and s-wave 
superconductivity^. The electronic specific heat was 
calculated 8 and quantitatively fit the experimental data 
quite well, as well as several other properties, including 
the nuclear relaxation ratei 

In this paper we show the results of a nuclear relaxation 
rate \jT\ calculation from a different starting modelifi 
in which there is p-wave pairing in the superconduct- 
ing channel. This result is then compared to the s-wave 
calculation and the experimental data to determine a sig- 
nature of the much disputed pairing symmetry. We find 
that there is a small signature around T c in the observed 
data of a small Hebel-Slichter peak that resembles the 



s-wave curve, but that this comparison is inconclusive to 
verify either type of pairing. 

A triplet state can be described by the order pa- 
rameter A a/3 (k) = (c k ^ a Ckjj) = [i{d(k).a)a y ] a f3, where 
a = xa x + ya y + za z denotes the usual Pauli matri- 
ces, and the basis of symmetric matrices iaa y is used 
to represent odd angular momentum pairing. The three- 
dimensional complex vector d(k) fully characterizes the 
triplet pairing state. For simplicity, the easy axis of mag- 
netization is assumed to be in the z-direction. Because 
of the pair-breaking effect of the strong exchange field, 
only the Cooper pairs with equal spin states will survive. 
In the case of equal spin pairing we can write vector d 
in the form d = (d x , d y ,0). Denoting A± = d x ± d y , the 
superconducting order parameter (SCOP) becomes^ 
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An attractive pair-forming interaction V is assumed, 
and the coexistent Hamiltonian is then 

Hfm+sc = ^2 £kC k,Q C k,a - I drs(r) • s(r) 
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The equations for the SCOP have the usual BCS-like 
form: 

A-(k) = zlVv(k l k') 1 - 2/ ^ k '»A-(k') (3) 
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where f(E) is the Fermi-Dirac distribution function, and 
A+ and A_ are the SCOP for the up-up and the down- 
down spin pairing states respectively. The equation for 
the magnetic order parameter (MOP) is 



M = 




2/(£_)] 4[1 - 2f(E + )} 
2E + (k) 



(5) 



where e k r,i = £k ± Eas. l|3l4l5|l are coupled with 

each other via the quasiparticle spectrum E±(k) = 

( ekT ^) 2 + |A ± (k)p. 

In order to illustrate the interplay between the MOP 
and the SCOP, we solve eas. l|3l4l5l) self-consistently for 
the simplest case of a spherical Fermi surface, assum- 
ing that the SC pairing strength is in the p-channcl, 

V(k,k') = ffi=i(fc,A:')Em=-i^(k)iT m (k'), and has 
the BCS-like form, i.e., g\{k, k') vanishes everywhere ex- 
cept the narrow region near the Fermi surface and does 
not depend on the exchange interaction /. We have used 
A±(k) = AoiY^ 1 . The SCOP will have point nodes at 
the poles of the Fermi surface. 

The order parameters, M and A, have dependencies 
such as M = M(g,I,T) and A = A(g,I,T) at T ^ 0. 
After solving these three coupled integral equations nu- 
merically we study the characteristic behavior of the or- 
der parameters. 

The variation of the transition temperature of the 
ferromagnetic and the superconducting transitions as a 
function of / has been studied previously 1 - , and here we 
mention the results. For / < 1.0, which corresponds 
to the paramagnetic phase M = 0, the superconducting 
transition temperature for the up-up, T c ,, and down- 
down, T c _ , spin paring is equal, implying that A + = A_. 
This corresponds to the planer (unitary) phase of the 
triplet pairing. But for I > 1.0, there is a finite mag- 
netization and the transition temperatures T c+ and T c _ 
will not be equal, implying that A + ^ A_. Upon fur- 
ther increase of /, and the subsequent increase of M, T c+ 
increases and T c _ decreases to zero. Therefore, the finite 
magnetization induces the non-unitary type of supercon- 
ducting triplet state. 

Here, we study the temperature dependence of M, A + , 
and A_ for fixed values of / and g. For / < 1.0, in the 
paramagnetic phase, the variation of A + and A_ is iden- 
tical to that of the mean field result A cx (T c — T) 2 as in 
the BCS case. Similarly, the temperature variation of M 
in the absence of superconductivity, i.e., A + = A_ = 0, 
follows the trend as is predicted by mean field theory 
(solid line in Fig.lJTJ). We show the variation of M, 
A + , and A_ as a function of T for the chosen values 
of I = 1.018 and g — 0.95 (coexistence of magnetism and 
both equal spin pairing triplet state superconductivity) in 
Fig. Q . The temperature axis has been normalized with 
the ferromagnetic transition (Curie) temperature T m to 
show the relative magnitudes of the corresponding tran- 
sition temperatures. We can tune values of I and g so 
that the transition temperatures will be comparable. 

We see in Fig-Q that the SCOP decrease with tem- 
perature. The MOP decreases as well for T > T c+ , 
but for T < T c+ we see some new features. If A + > 
A_ ^ 0, for T < T c+ , M (dashed line) is always big- 
ger than its value in the pure ferromagnetic state when 
A + = A_ = (solid line). When there is a finite A + 
and A_, M increases with the increase in T from to 
T c _ , and then decreases sharply in the range of tempera- 




FIG. 1: Order parameters in this model as a function of T. M, 
A+, and A_ are labeled accordingly. For M, the solid curve 
refers to the pure ferromagnetic state, A + = A_ = 0, the 
dashed line refers to M when A+ > A_ 7^ 0, and the dotted 
line when A_ = and A + 7^ 0. The y-axis is a unitless 
measure and the x-axis is scaled by the Curie temperature 

(Tm). 



ture T c _ < T < T c+ (dashed line). The increase in M for 
< T < T c _ is due to the decreasing A_ as the temper- 
ature increases. On the other hand, the sharp decrease 
in M for T c _ < T < T c+ is not only due to the temper- 
ature, but more strongly to the decrease of A + . Again 
we see that for < T < T c _ , M is bigger if A_ = 
(dotted line) than when A_ ^ (dashed line). These 
results are in agreement with the general intuition that 
the magnetization has a cooperative relation with the up 
spin pairing state, and a competitive relation with the 
down spin pairing. Our ultimate goal is to study the 
variation of the relaxation rate 1/Ti with temperature 
for this p-wave state, so next we derive the equation for 
the density of states as a function of excitation energy. 

To derive the expressions for the density of states we 
use the usual relation, N%(E) = J dpS(E ~ E°), 

where a refers to the up spin or down spin fermions. 
Solving in the standard way using the property of the 
5— function, we get, 



N^tE) = At / dOsind (A+ + A~) , 
8tt 2 J q v ; 
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FIG. 2: Density of states as a function of the excitation energy 
for each spin species. The dotted line corresponds to the 
normal state at the Fermi level. Inset A corresponds to M = 
0, B to M ^ for which down spin pairing is finite and C to 
M^O for which down spin pairing is completely suppressed. 
Note that the density of states of the up and the down spin 
fermions are not identical when M/0. 



The + sign is for the | spin fermions and the — sign 
is for I spin fermions on A±, respectively. These two 
expressions converge to the density of states of the nor- 
mal state at the Fermi level, N(Q) = m % f , in the 
limit of all OP equal to zero. For M = 0, the ex- 
pressions result in the density of states of the equal 
spin pairing triplet state superconductor, N(E)/N(Q) = 
E/2fd6sm9/ v / (E 2 - A 2 ). The important effect of the 
internal magnetization M is to modify the value of the 
density of states. For increasing M , we observe an in- 
crease in Nj (E) and a decrease in Nj (E) . 

We plot the T — density of states (DOS) scaled to 
the normal state value for positive excitation energies in 
Fig. © • The inset A represents the density of states when 
M = 0, which is just a superconducting DOS and both 
spin components fall on the same line. The divergence 
in the density of states is at E = A + = A_. The main 
body B of Fig.© shows the effect of the magnetization 
on the density of states when the magnitude of M is such 
that A_ 7^ 0. We see that due to the finite magnetiza- 
tion, the density of states curve now splits for the up and 
the down spin fermions. The diverging density of states 
at the lower energy is at E = A_ and that at the higher- 
value is at E = A + . For the value of I at which M is 
large enough to completely suppress the down spin pair- 
ing, the density of states for the down spin fermion will 
be constant as is shown in inset C. These modifications 
in the density of states due to the presence of the internal 
magnetization will manifest themselves in the nuclear re- 
laxation rate of the system. 

In general \/T\ in the superconducting state is related 
to the (DOS) in the superconducting state, N S (E) as 



followsii: 

1 r°° 

-ex 2 / (N s (E)ff(E)(l-f(E))dE, (8) 

where f(E) is the Fermi-distribution function. In the 
presence of M the functions will not be the same for the 
up and the down spin fermions as mentioned above, so 
we split the integral for different spin species. Hence, 

1 r°° 

-ex / (N±(E)ff(E-)(l-f(E-))dE 

Ti J ° o) 

+ / (Nj(E)ff(E + )(l-f(E + ))dE. 
Jo 

First, we choose the interaction parameters I and g 
such that there is both the up-up spin and the down-down 
spin pairing (the corresponding transition temperatures 
are not equal) and study the variation of \jT{T with re- 
spect to temperature. The results have been presented 
in Fig.©, where the temperature has been normalized 
with the transition temperature for the up-up spin pair- 
ing. We choose I and g such that T c+ /T c _ ~ 2.5. In 
Fig.© we see that l/T x T decreases for T < T c+ . This 
decrement characterizes the onset of the spin triplet pair- 
ing. With the decrease in temperature the rate contin- 
ues to decrease, but not as strongly as it would for the 
M = shown by the dotted line. The dotted 

line is just the p-wave superconducting case, and when 
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FIG. 3: 1/T\T vs. normalized temperature with respect to 
the transition temperature of the up spin Cooper pairs, T c +. 
The dotted line is the p-wave superconducting case, with no 
magnetization. The top solid line is for the case when there 
is no pairing of down spins, where the only transition is at 
T/T c + = 1 when the system becomes an up spin paired su- 
perconductor, and saturates at T = 0. The bottom solid 
curve shows a second transition when the down spin Cooper 
pairs are formed and another gap A_ opens up, and is fully 
suppressed to zero. The values of I and g are chosen so that 
T c+ /T c . = 2.5. 
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FIG. 4: Comparison of theoretical and experimental data. 
The temperature axis is scaled by the SC transition temper- 
ature. The experimental data— was taken at 13kbar: at this 
pressure there is an observed superconducting transition at 
T ~ 0.5K 

the gap opens up, the rate is suppressed quickly. 

On the other hand, when we have a finite magnetiza- 
tion and only superconducting pairing of the up spins, we 
refer to the top solid curve of Fig.©. In this case, the 
transition at T — T c + results in a decrease of the relax- 
ation rate due to the presence of the gap A + , but the rate 
does not decrease as sharply as the unpolarized case be- 
cause there is a finite magnetization. Furthermore, there 
is a normal fermionic component in this case since the 
down spins have not paired (due to our choice of the in- 
teraction parameters) , and thus this finite density at the 
Fermi level allows for a nonzero relaxation rate at zero 
temperature. This follows from the fact that the entire 
Fermi surface is not gapped. 

Finally, the bottom solid curve of Fig.© is calculated 
when we allow for the down spins to pair. This results in 
a second superconducting transition at T c _ and the open- 
ing up of a second gap A_. The relaxation rate goes to 
zero when all the states are frozen out at zero temper- 
ature as expected. Note that there is no Hebel-Slichter 
peak for the p-wave case. 

We now compare these results with the s-wave singlet 
nuclear relaxation rate and the experimental data. We 
have received data from the Kotegawa groupi 2 on XJG&2- 
The calculation of the singlet has been carried out in a 
prior publications., where the details can be found. In 
Fig.© we calculate 1/Ti versus temperature and show 



all three curves. The experimental data was taken at 
Yikbar. We emphasize that at this pressure and tem- 
peratures, there is coexistence of superconductivity and 
magnetism. Also, the theoretical triplet curve is for the 
case of one single transition at T = T c corresponding 
to the up spins forming Cooper pairs. The down spins 
do not pair in this case and behave basically as normal 
fermions. 

Several features are seen in Fig.©. First of all, the 
data is not very consistent near the critical temperature, 
due to the difficulties in the measurement. However, a 
small peak below T c is noticeable, which follows the trend 
of the s-wave case of the solid line. This very well could 
be the reduced Hebel-Slichter peak that was predicted 
due to the presence of the ferromagnetism which pro- 
duces a finite density of states at the Fermi level due to 
the existence of gapless fermions^. 

This p-wave curve is calculated without down spin 
pairing, and referring to Fig.© where \jT\T is constant 
at small temperatures, 1/Ti goes to zero with T as ex- 
pected. The low temperature behavior of the data seems 
to agree well within the triplet scenario. We do not ob- 
serve in the experimental data a second transition and 
a second gap opening up from the pairing of the down 
spins. From our theoretical analysis, we would expect 
this to happen at a temperature of around T = 0.018X, 
however the experiments have only been carried out down 
to T = 0.04^. 

In conclusion, we have studied a p-wave model of su- 
perconducting pairing in the background of itinerant fer- 
romagnetism and calculated the nuclear relaxation rate 
1/Ti. Since there is still much debate as to the type of 
pairing in these systems, we feel a direct comparison be- 
tween what is expected for different types of pairing and 
experimental data is the most important goal. We see 
that around the critical temperature, the s-wave model 
agrees well with the experiment, as it does in the specific 
heal£, but at lower temperatures there are discrepan- 
cies with the theory. The model does make clear that 
the magnetization suppresses strongly the superconduct- 
ing signatures, which is certainly being observed experi- 
mentally, and a closer look at other electronic properties 
in conjunction with careful experiments is what is need 
to determine the nature of the superconducting pairing 
symmetry. 

We thank the Kotegawa group for providing an early 
draft of their paper— on UGe^- This work was done with 
the support of DOE/DEFG0297ER45636. 
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